Island formation and dynamics of gold clusters on amorphous carbon films 
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Samples of Au clusters deposited by laser ablation on an amorphous-carbon substrate are in- 
vestigated. After a few months storage at room temperature the initially statistically distributed 
clusters are found to be collected in agglomerates consisting of larger clusters embedded in an Au 
film typically covering areas of size 25x70 nm^. The Au film is determined to be probably 4 to 
8 monolayers but at most 7 nm thick. Evidence is found that a number of clusters consisting of 
less then 50 atoms are pinned at intrusions of the substrate. These results were derived using high 
resolution transmission electron microscopy and off-axis holography measurements to characterize 
the agglomerates as well as the substrate. Monte-Carlo simulations were performed to model the 
film formation process. To this end the substrate-Au interaction was determined using density 
functional calculations (GGA) while the Au-Au interaction was modeled with effective many body 
Gupta potentials. The film formation can be understood as diffusion and fusion of clusters of in- 
termediate (50 < N < 300 atoms) size. Larger clusters are more stable at room temperature and 
remain adsorbed on the Au film. 



I. INTRODUCTION 

I The dynamic behavior of clusters deposited on a sub- 
strate has been studied for several decades ii2i2i^iS4& A 
continuously growing number of potential applications 
in electronics^'- and catalysi a^i^° demands defined arrays 
of nanoparticles. However, such arrays tend to mini- 
mize their energy by minimizing the total surface via 
the growth of larger particles at the expense of smaller 
oncs.^ An example for such a process is Ostwald ripen- 
ing, where the larger cohesive energy in larger particles 
leads to matter transport away from small particles in the 
presence of a finite partial pressure of the constituent in 
the environment surrounding the clusters>ii*i^ Therefore, 
investigations concerning the stabilit}ii2*i^ of deposited 
nanoparticles are of considerable interest with respect to 

[ potential nano- technological applications. 

Moreover, the thermodynamioi^ii^ and chemicaUi 

|. properties of small metal clusters themselves have en- 
joyed a large interest over the past years. Their prop- 
erties differ from those of the bulk material raising the 
fundamental question about the statistical mechanics of 
finite systemsii^ Since the presence of a substrate al- 
ters these propertiesSii^ detailed investigations of the 
substrate-adsorbate interaction's are required. 

In this paper we present results from transmission 
electron microscopy (TEM) investigationsSiiS^ concern- 
ing the time- and temperature-dependent behavior of 
Au clusters deposited by laser ablation technique on 

' amorphous-carbon (a-C) films. The reasons to choose 
an amorphous substrate are threefold. Firstly, the non- 
crystalline structure of the substrate allows for the visu- 
alization of the crystal structure of the adsorbed clusters 
in high resolution TEM (HRTEM) experiments. Sec- 
ondly, amorphous carbon is mechanically much more 
stable than highly oriented pyrolytic graphite (HOPG) 
and consequently much more suited for technical appli- 
cations. Thirdly, as opposed to the conductor HOPG, 
amorphous carbon is semi-conductingS2iSi and as such is 
expected to be less influential on the electronic proper- 



ties of the adsorbate. Furthermore, our observations were 
made on samples that were stored under normal atmo- 
sphere for several months enhancing their relevance for 
possible technical applications. 

The central observation reported in this paper is the 
formation of islands on a time scale of several months 
after the deposition of the Au adsorbate. Typical TEM 
images of the samples are shown in Fig. ^ Panel (a) is a 
sample one day after its preparation. The inset shows the 
size distribution of the clusters. Clusters smaller than 1 
nm in diameter cannot be resolved unambiguously be- 
cause of their low contrast on the amorphous substrate. 
Panel (b) shows a typical island formed after the sam- 
ple has been aged for four months. The reproducibility 
of the results is assured since eight samples of different 
coverage prepared on different days stored in individual 
sealed containers show the same phenomenon. 

While the formation of small islands on short time 
scales has been observed previously for larger size- 
selected clusters for Co on microgrid substrates,^ for Au 
on amorphous carbon as well as for Ag on graphite,^ 
the new features here are dark areas typically covering 
25x70 nm^ underlying the larger clusters. As described 
in this paper in detail we were able to identify these ar- 
eas as islands consisting of Au films of roughly 4 to 8 
monolayers thickness with immersed larger Au clusters. 
The Au films show areas of crystalline structure with es- 
sentially two different lattice plane orientations as shown 
in the circled areas in Fig. ^c) . A close analysis of the 
properties of the Au films including a holographic deter- 
mination of its thickness is given is Sec. llVI We arrive at 
the conclusion that the islands are formed by coalescing 
small clusters with diameters of < 2 nm. 

Since the observed island formation was unexpected 
on the amorphous substrate, no monitoring of the dy- 
namics of the formation of the islands — as desirable — was 
done. Awaiting the preparation of new samples, which 
are needed for the observation of the formation process, 
the present paper aims to give a state-of-the-art analy- 
sis of the system at hand. Based on these findings the 
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FIG. 1: TEM images of Au clusters deposited on an a-C 
substrate by laser ablation (3000 shots): (a) sample one day 
after its preparation. The inset shows the histogram of the 
size-distribution of the clusters. Clusters smaller than Inm in 
diameter cannot be resolved unambiguously within the TEM 
images, (b) Typical island of higher contrast formed on a 
sample four months after its preparation, (c) Detail of the 
island displayed in panel (b). The areas circled show regions 
exhibiting crystalline structures with differently spaced lattice 
planes. 



time-analysis will be performed once new samples are 
available. 



Outline of the paper 

In Sec. m we describe the experimental [sample prepa- 
ration and TEM] as well as theoretical [density functional 
theory (DFT) and Monte-Carlo (MC) simulations] meth- 
ods applied. 

Since the main trade-off of the amorphous substrate is 
its less well defined surface we devote Sec. lIIII to develop 
a model of the carbon film that is consistent with the 
experimental observations. To this end holographic im- 
ages of the substrate film are analyzed fSec. lIII First 
principle DFT calculations allow for the determination 
of the substrate- Au interaction fSec lIIIB)) . 

In Sec. llVl the observed islands are analyzed in detail. 
Their thickness is determined holographically fSec lIVB)) 
and their stability tested in a heating experiment (Sec. 

In Sec. the formation of the islands and the shape 
of the Au clusters on their surface is modeled with MC 
simulations followed by the conclusion Sec. I VII 

II. METHODS 

A. Experimental: substrate and sample 
preparation 

The commercial a-C substrate films were produced by 
evaporation in a carbon arc by Arizona Carbon Foils 
and distributed by Piano GmbH as type S160. The 
films arc mounted on a 200 nm mesh Cu grid for sup- 
port. The film thickness is given by the manufacturer 
as dsubs = 10 — 12.5 nm with a density of Psubs ~ 2.0 
g/cm^ corresponding to a mass of 2 — 2.5yLtg/cm^. The 
density of the substrate is closer to that of graphite 
(Pgraph = 2.267 g/crn^) than to that of diamond (pdia = 
3.515 g/cnr^) suggesting a structure that contains re- 
gions with trivalent coordination^ The similarities in 
the electronic structure between a-C and graphite are 
supported by electron energy loss spectroscopy (EELS) 
measurements )S4 which show spectra that differ signifi- 
cantly from those of diamond. Amorphous-carbon films 
produced with similar methods and of similar density 
have been found to be semi-conductingi22i2^ 

Note that, since psubs < Pgraph and since in the amor- 
phous material some tetravalent C with higher local den- 
sity must be present, there must exist areas with very low 
density or even voids for compensation. This assumption 
is supported by the results presented in Sec. IIIII 

The Au clusters were collected from the primary beam 
of a laser vaporization cluster source, which has been de- 
scribed elsewhere in detailiSSiSi In brief, the laser vapor- 
ization cluster source is a variant of the Smalley-deHeer- 
^ypQ28j29 gg(;^p optimized by Heiz'^" for high yield. The 
source is equipped with a rotating gold disc target with 
a diameter of 50 mm which is sealed with a Tefion gasket 
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against the source block. A pulsed laser (Neodym-YAG, 
Continuum, 532 nm, 30 Hz repetition rate) is focused 
through a nozzle onto the target. A pulsed valve (General 
Valve, 5 bar backing pressure of He) which is synchro- 
nized with the laser quenches the evaporated atoms into 
clusters which expand through the nozzle and a skimmer 
into an oil diffusion pumped vacuum chamber at 10^^ 
mbar. Au clusters (and atoms) are deposited without 
further mass selection onto a TEM grid placed in the 
primary beam in a distance of about 40 cm from the 
nozzle. 

The size distribution of the Au particles is shown in 
the inset of Fig.ma). Since the only directed acceleration 
of the clusters is the expansion into the vacuum of the 
cluster source the impact energy can be considered small 
enough to avoid intercalation of adsorbate and substrate. 

The investigations were conducted on a batch of 8 sam- 
ples prepared on different days. Four samples were pre- 
pared with 300 ablation shots, four with 3000 ablation 
shots. They were stored individually under air in sealed 
containers. Experiments were carried out to assure the 
absence of contamination on the samples fSec. II V|). All 
samples show the same island formation assuring the re- 
producibility of the results. 

B. Experimental: transmission electron microscopy 

The TEM was carried out with a Philips CM200 
FEG/ST electron microscope at an energy of 200 keV 
equipped with a Noran Ge detector system for energy- 
dispersive X-ray analysis (EDX). A Gatan 652 Double 
Tilt Heating Holder operated by a Gatan 901 SmartSet 
Hot Stage Controller was used to perform the in situ an- 
nealing experiments. 

Transmission electron holography was carried out us- 
ing a MoUenstedt biprism installed in the selected-area 
aperture holder of the microscope. The electrostatic po- 
tential of the biprism wire was close to 150 V. The im- 
ages recorded were analyzed using the phase shift of the 
(OOO)-beam of the first hologram sideband. Data analysis 
was performed using the DALI program package;^ which 
was extended for the reconstruction of holograms. The 
details of the reconstruction of holograms are outlined by 
Lichte and LehmanniS^ 

Figure |21 shows the three-dimensional visualization of 
the phase shiftH with respect to the vacuum A0 observed 
on a 28.2x28.2 nm^ surface segment of a typical sample. 
In the far right corner a part of an island as shown in 
Fig. n^b) is visible. In the front left corner the linear 
decrease of the thickness of the substrate near the sub- 
strate edge is visible. All holographic images were taken 
at substrate edges in order include a section of the vac- 
uum for calibrationiSiiS^ A closer analysis of the data is 
given in Sec. IIII Al concerning the substrate and in Sec. 
IIV Bl concerning the Au adsorbate. 

The relation between the change in phase shift A(p ob- 
served in the holographic images and the corresponding 
thickness d of the observed object can be derived from 
scattering theory^! and is expressed via the standard for- 




FIG. 2: Holographic phase shift A(j) of a segment of the a- 
C substrate with part of an island as shown in Fig. ^b) in 
the far right corner. The image has the size of 28.2x28.2 
nm^. The z-axis value is proportional to the observed phase 
shift. In the foreground the almost linear decrease of the a- 
C film thickness towards its edge is visible. (For alternative 
representations of the data see Figs.\^and\^) 

mula 



where the interaction constant depends on the acceler- 
ation voltage resulting in Ce — 7.29 x 10^ rad/(V m) 
in our case and where Vq denotes the effective internal 
potential of the sample. 

For bulk gold the theoretical values^SiM for Vq are 
VAu.thoo = 28±2 V, while experimental values^^ are given 
as Vau.cxp = 22± 1 V. For small Au clusters on Ti02 sub- 
strates a sharp increase of Vq for clusters smaller than 4 
nm was observed with values of up to V^u. cluster ^ 

50 

V for clusters smaller than 2 nm>i2i We confirm these 
findings in Sec. IIV Bl 

For the a-C substrate films the situation is less consis- 
tent. We observe phase shifts of the substrate film with 
respect to the vacuum of ~ 1.9 — 2.8 radians. Associating 
these phase shifts with the numbers of the manufacturer 
for the film thickness of 9.5 to 14 nm we obtain a poten- 
tial of Fsubs = 27.4 V. 

On the other hand, the experimental value for a-C films 
of the same density as ours has been given by Harscher 
and Lichte'^^ as Vhl = 10.7 V. This number compares 
satisfactorily with the values for graphite summarized by 
Sanchez and Ochandc'^'' of Igraph.thco ~ 12 — 15 V (the- 
ory) and Fgraph.cxp ~ 11 — 13 V (experiment) since the 
ratio of the potentials ^HL/Vgraph ~ Psubs/Pgraph corre- 
sponds roughly to that of the densities, i.e., the inner po- 
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tential scales roughly with the density of the material. A 
similar argument holds for the values for diamond, where 
Vdia,thco 16 - 23 V (theory) and Vdia^cxp 15 - 21 V 
(experiment). 

The obvious discrepancy between the values of Ksubs 
and Vhl cannot be resolved here. Possible origins of the 
discrepancies are (i) the underestimation of the substrate 
thickness by the manufactures^^i^ and (ii) an increase of 
the internal potential for thin films similar to the increase 
of the internal potential in small Au clustersii^i Note that 
the latter effects would not have been observed for the 
thicker films with d > 35 nm in the investigations by 
Harscher and Lichtei^ We will discuss the implications of 
the discrepancy for the present work where appropriate. 

C. Theory: density functional theory calculations 

We have calculated the interaction between Au atoms 
and the carbon substrate using first principles DFT 
calculations in the generalized gradient approximation 
(GGA).""^ For all calculations we used the highly accurate 
Projector Augmented Wave method^ as implemented in 
the VASP electronic structure programi^ 

For the carbon substrate we used both a graphite sur- 
face as well as a model a-C surface consisting of 150 
carbon atoms4^ The Au-substrate interaction was calcu- 
lated for a number of different positions of an Au atom 
relative to the substrate surface. For each position the 
Au-substrate interaction was mapped out by holding a 
Au atom fixed at different heights above the surface while 
the carbon atoms at and near the surface were allowed to 
fully relax. As a reference point a Au atom at a distance 
of 9.5 A from the surface was used fSec. IIII B|| . The re- 
sulting Au-substrate potentials were averaged and fitted 
to a modified Lennard-Jones potential and used as input 
in MC simulations described in the next section. 



D. Theory: Monte-Carlo (MC) simulations 

In order to develop a microscopic understanding of the 
dynamics of the Au clusters on the substrate we simulate 
the system with the canonical Monte-Carlo method. A 
standard Metropolis algorithm is employed^^'^^'^ with 
an update after each random displacement of an Au atom 
within an interval [0, dmax] in all spatial dimensions, dmax 
is set to yield an MC acceptance rate of 50 to 60 %. 
The resulting temperature dependence is roughly dmax oc 
\/T . The boundary conditions are imposed by a hard wall 
cube with linear dimension L^^ Ly, and L^. 

The Au-Au interaction is modeled via the many-body 
Gupta potentiali^ (GP): 

N N 

^ I 

- T.JJ2^' e-^.C^-./^o-i) . (2) 



The distances = \ri — rj \ are measured in units of the 
bulk first-neighbor distance ro = 2.885 A. The indices 
i,j G {1, . . . , N} label the Au atoms at positions and 
Tj, respectively. The parameterizations for Au as found 
in the literature^ i.e., A = 0.2061 eV, ^ = 1.790 eV, 
p = 10.229, and q = 4.036, has been determined to match 
the bulk elastic constants and the surface contraction. 

The Au-substrate interaction is modeled by two- 
particle interactions in the form of generalized Lennard- 
Jones 6-12 and 3-6 potentials. Their derivation is dis- 
cussed in detail in Sec. IIII Bl 

Runs on the atomic level have been performed for up 
to 10* updates per atom for up to 600 atoms. Runs 
have been stopped when a metastable configuration is 
obtained that does not evolve anymore on a reasonably 
accessible time scale. Depending on the size of the sys- 
tem metastable configurations are attained after a few 
minutes {N = 55) or a couple of days (Au double layer 
formation with embedded clusters with N — 600). Runs 
where performed up to three weeks to assure that the life- 
time of the metastable configuration is at least an order 
of magnitude larger than its formation time. Run times 
for particular cases are given with the resutls in Sec. |3 

In order to simulate larger systems with many clusters 
we derive in Sec. IV Al an effective cluster-cluster poten- 
tial. Runs are then performed with a few hundred clus- 
ters to illustrate the cluster diffusion and fusion process 
in Scc.rvTH 

III. SUBSTRATE CHARACTERIZATION 

The observed formation of the islands as shown in 
Fig. n implies that the substrate surface is free of strong 
pinning centers. The same conclusion can be drawn from 
oval shapes of the islands with smooth boundaries. More- 
over, the pattern formation of the Au film discussed in 
Sec. IIV Al [Fig. ^c)] strongly suggest the presence of 
smooth areas on the substrate surface. 

Since this observation is not ad hoc intuitive for an 
amorphous substrate we investigated the substrate with 
off-axis holography experiments. In order to obtain 
a qualitative understanding of the underlying mecha- 
nisms of the island formation we determine the substrate- 
adsorbate interactions via DFT calculations. The latter 
are used as input parameters for MC simulations which 
reproduce the experimental observations. 

A. Holography 

When imaging vacuum with TEM off-axis holography 
the phase shift observed is slightly fluctuating due to im- 
perfections of the biprism wire, aperture, and noise in- 
duced by the CCD camera. We analyzed these vacuum 
fluctuations in order to distinguish them from the sample 
signal. The inset in Fig. |SIa) shows a typical image of 
the vacuum phase shift (/)vac after subtraction of its mean 
value. The main graph of Fig. OJa) shows the normal- 
ized distribution w^g^d^^'vac) of the vacuum fluctuations 
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nm. Such narrow intrusions with a depth of up to 1/5 
of the substrate thickness would make the substrate very 
unstable. Moreover, the edges of the intrusions would in- 
evitably lead to strong pinning centers similar to those of 
nanopits in HOPG,''^'^" which are inconsistent with the 
observed island formation. 

Since the intrusions are only 2 to 4 C-C bond lengths 
wide, it appears likely that they are capped in the pro- 
duction process thus stabilizing the structure. The in- 
trusions are then left as voids in the film. The effective 
potential is strongly reduced at the position of the voids 
leading to the observed phase modulation. At the same 
time the cap yields a rather smooth surface consistent 
with the observed island formation. The C atoms sur- 
rounding the voids have less nearest neighbor atomS)^ 
which is consistent with the measured EELS spectr£k24 
indicating largely trivalent coordination. In the absence 
of methods for a more precise determination of the lo- 
cal structure of the substrate we use this scenario as a 
working hypothesis. 

Finally, the histogram in Fig. shows the normal- 
ized phase shift distribution of the a-C substrate film with 
the Au adsorbate, albeit at positions without any of the 



observed islands. (The latter are discussed in Sec. lIVB|l . 
The Gaussian fit from Eq.|3|as shown by the full line has 
a width of CTgubs — 0.164 rad which is 25% larger than 
(Ta_c- The difference is readily interpreted as induced 
by small Au clusters pinned at small intrusions of the 
substrate. Since the pinned clusters are smaller than the 
TEM resolution limit their diameters must be smaller 1 
nm or < 50 atomsi^ 



FIG. 3: (Color online) Phase shift distributions as a measure 
of fluctuations of the phase. Panel (a) shows the result for 
the vacuum, panel (b) shows the substrate Elm without ad- 
sorbate, and panel (c) shows the distribution for the substrate 
film with Au adsorbate albeit at positions without islands. 
Please note the different scales of the x axes. Histograms 
are obtained by binning and normalization of the shifted raw 
data, full lines are Gaussian fits to Eq. The insets show 
samples of the measured phase shifts. The data used to ob- 
tain the histograms amounts 5 to 10 times the date shown in 
the insets. 
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from which we obtain a width of the vacuum fluctuations 
of (7 = CTvac = 0.052 rad. 

The histogram in Fig.|3^b) shows the normalized phase 
shift distribution of the a-C substrate film without Au ad- 
sorbate. The Gaussian fit from Eq.|3|as shown by the full 
line has a width of tTa-c = 0.13 rad. Since da-c > 2crvac 
the fluctuations must reflect properties of the substrate. 
Using Eq. ^ and the value of T^ubs — 27.4 V the change 
in thickness corresponds to Adgubs = 2(Tsubs/(C'E^ubs) = 
1.3 nm, for Vhl = 10.7 V even Adgubs = 3.3 nm. The 
graph in the inset Fig. I^Jb) reveals that phase fluctua- 
tions of comparable depth occur on length scales of 0.3 



B. Surface adsorbate interaction: DFT calculations 

A quantitative description of Au/a-C substrate inter- 
action is needed for the MC simulations described in Sec. 
El To this end we studied the interaction of individual 
Au atoms with graphite and a model a-C surface using 
ab-initio DFT calculations. 

For the calculation of the Au/graphite-surface interac- 
tion we used a supercell consisting of 2x2 graphite unit 
cells in the x-y plane and four carbon layers in the z- 
direction. A large vaccuum distance of 19 A was usedi^ 
Calculations were performed for Au in the top, bridge, 
and hollow sites. Using as reference a Au atom at the 
center of the vacuum we calculated binding energies of 
0.065 eV, 0.062 eV and 0.041 eV, respectively The Au- 
graphite potential for Au in the top-site as a function of 
Au atom surface distance is shown by line "g" in Fig. 
2] The result of the DFT calculations was fitted with a 
modified Lennard- Jones potential 



^J, graph — Vlj,0 
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1.11 
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(4) 



where Vlj,o is the binding energy and ro = 2.885 A. 

For a more realistic description of the a-C substrate 
we have constructed a model a-C surface. To ensure 
that sp^ bonding will be dominant^'* and to obtain a 
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FIG. 4: (Color online) Au/graphite (squares) and Au/a-C 
substrate potentials (circles) as calculated within the DFT 
approach. The inset shows the corrugated surface layer of 
the a-C model with the locations of the calculated potentials 
(a)-(d) indicated by arrows. The fitted potentials used in the 
MC simulations are shown as solid lines for graphite (g) and 
the a-C model substrate. Dashed lines are guides to the eye. 



Au atom to C sites lying in such a depression are consid- 
erably larger as compared to the graphite surface (cases 
b,d in Fig. 0)). An Au atom in the center of the depres- 
sion can bind to several C atoms and the binding energy 
is smaller but of the same order as the Au-Au binding 
energy (case c in Fig. 0}. 

To model the binding energy of Au atoms in clusters 
with a contact area larger than the surface intrusions 
the resulting Au/a-C substrate potentials were averaged 
using approximate relative surface areas as weights and 
fitted to a modified Lennard- Jones potential 

y ( 297 34.5 \ 

KLJ,n.ea„ - 1/L,I,0 ^^^^^ _ ^ 2)12 (z/ro - 1.2)6 ) ' 

(5) 

where Vlj,o is the binding energy and tq = 2.885 A. Using 
weights of c:d:b:a=l:4:9:16 results in an average binding 
energy of Vlj.o = 0.34eV. The MC simulations in Sec. Ivl 
do not depend qualitatively and only little quantitatively 
on the specific value in a range of 0.2 < Vlj.o/cV < 0.4 
and hence an average value of T^j^o = 0.3eV was adopted. 

The value of Flj.o ^ 0.3 eV is only appropriate for 
sufficiently large clusters. Small clusters with contact 
areas smaller than the intrusion size have larger weight 
of the strongly binding intrusion center c and are easily 
pinned consistent with the observations discussed in Sec. 

luSI 



simple model consistent with the observed Au-film for- 
mation we use graphite as our starting point. We use a 
supercell consisting of 4x4 graphite unit cells in the x-y 
plane and four carbon layers in the z direction. Addi- 
tional carbon atoms were added in localized regions but 
otherwise randomly inbetween the graphite planes so as 
to create voids consistent with the observed holographic 
phase fluctuationsi^ 

The resulting supercell of 150 carbon atoms was re- 
laxed into a local minimum. Two main effects of the 
relaxation can be observed: The carbon structure ex- 
panded strongly in the z-direction and the graphite lay- 
ers buckled due to the uneven distribution of interstitial 
carbon atoms. The density of the final a-C model struc- 
ture is p ^ 2.0 g/cm'^ in good agreement with the density 
of the a-C substrate as given by the manufacturer. The 
surface layer of the final structure is smooth but has a 
corrugation of d = 0.9 A over a distance of a few C-C 
bondlengths (Fig.^I). 

We have calculated the Au/a-C substrate potential at 
several positions roughly along a line through a depres- 
sion in the surface. Since for the graphite surface the 
Au on top position is the preferred bonding site, only Au 
on top sites were considered in the a-C case. High lying 
C surface atoms retain their graphite character and are 
predominantly 3-fold coordinated. Accordingly we find 
the binding energy of an Au atom is still comparable 
to the binding energy on graphite (case a in Fig. C 
atoms lying in a depression of the surface have additional 
bonds to the nearby C atoms below and are mostly four- 
fold coordinated. The calculated binding energies of an 



IV. LONG TIME-SCALE DYNAMICS 

The long time-scale rearrangement of heavy atoms on 
a-C substrates has been observed previously. ^'^ Here we 
present a detailed investigation of pattern formation of 
Au on a-C substrates. Figure ^^a) shows a typical TEM 
image of a sample prepared by deposition of 3000 shots 
of Au clusters on a 10 nm a-C film. The image displayed 
was recorded one day after the preparation of the sample. 
As expected, one observes a statistical particle distribu- 
tion. The histogram included shows the size distribution. 
Au clusters smaller than 1 nm in diameter could not be 
discerned from the 10 nm carbon substrate background 
owing to their low contrast. However, ion mobility mea- 
surements have revealed^ that clusters in this range of 
sizes are widely present in cluster ion beams. This sug- 
gests that a significant number of small clusters is also 
present on the substrate, albeit unresolved in the TEM. 

A. Island formation 

After the sample is kept at room temperature and in 
absence of inert conditions for four months, TEM stud- 
ies yield images as displayed in Fig. ^b). The previ- 
ously statistically distributed clusters are now collected 
in agglomerates referred to as islands. The vast major- 
ity of these islands is of similar size with typical areas of 
A ~ 25 X 70 nm^. 

In all islands an underlying area of higher contrast is 
observed, which is bordered by the outer clusters of the 
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islands. Using EDX-detection no elements except for Au 
and C are found in these islands. The binary alloy phase 
diagram of Au and C excludes the formation of Au-C 
alloys under the given conditions^'* and we conclude that 
the C signature stems from the substrate. The analysis 
given m Sees. HVBl HVTI and HVDl consistently suggest 
that the islands are formed of 4 to 8 monolayers of Au 
on top of the substrate. This conclusion is supported by 
the MC simulations of the formation process of such an 
Au film in Sec. |V| 

HRTEM images of the islands as depicted in Fig. ^c) 
show regions, in which crystalline patterns occur with lat- 
tice spacings determined as 501 pm and 310 pm within 
an estimated error of ±5%. The angles between these 
two spacings differ up to 6° from a rectangular configu- 
ration. It cannot be excluded that much larger areas of 
the Au films are crystallized forming such patterns be- 
cause a slight corrugation of the substrate may obscure 
its detection with HRTEM. 

While the exact determination of the origin of the 
aforementioned patterns from the present data is not 
possible, it can be speculated that the Au film ex- 
hibits a similar reconstruction as has been observed 
for Au(lll) surfaces by scanning tunneling microscopy 
(STM) measurements f^^i^ where superstructures with a 
similar modulation^ (y^ x lattice constant « 500 pm) 
have been observed. The presence of Au(lll) surfaces is 
also consistent with the MC simulations (Sec.0). More- 
over, reconstruction pattern have been observed for small 
islands of Au monolayers on HOPG in ultra high vacuum 
(UHV) by STM measurements but a direct compari- 
son is difficult because of the larger film thickness in our 
samples fSec. llV Bj) . 

We conclude that the higher contrast areas underlying 
the islands consist of a few monolayers of Au that form 
crystalline structures on at least parts of the substrate 
surface, which in turn must be sufficiently smooth. 

B. Au-film characterization 

In order to get information about the thickness of the 
observed films, we performed off- axis holography exper- 
iments. A three dimensional visualization of a segment 
of the substrate including a part of an Au island has 
been shown in Fig. [21 in Sec. Ill Bl Figure shows the 
same segment as an intensity plot, where lighter shades 
of gray correspond to larger phase shifts. The thick lines 
correspond to different scans investigated in detail. The 
circled areas indicate the positions of clusters adsorbed 
on the Au films that have been studied closely. 

The results for the three scans are qualitatively equiva- 
lent and for simplicity we focus the presentation on scan 
1. Figure |SJa) shows the corresponding phase profile. 
Four regions are clearly distinguishable from right to left: 
the rise of the substrate edge, the plateau of the sub- 
strate film, the rise of the Au-film edge, and finally the 
Au-film. The dashed lines show the profile after averag- 
ing the phase fluctuations. The vacuum level was deter- 
mined in a larger area a bit further away from the sample 




FIG. 5: (Color online) Holographic phase reconstruction im- 
age from the same segment as shown in Fig. as an intensity 
plot. Lighter shades of gray correspond to larger phase shifts 
with respect to the vacuum. Lines show scans that underwent 
closer investigation, circles indicate the positions of selected 
adsorbed clusters. 

(not shown). For the substrate edge the dashed line is a 
quadratic fit, for the Au-film edge it is a linear fit. The 
substrate plateau is given by a horizontal line at the av- 
erage value of phase in that region. Finally, the Au-film 
base line is obtained by adapting the average value of 
the fluctuation phase without the adsorbed clusters as 
outlined closer below. 

Panel (b) of Fig. shows the fluctuation of the data 
after subtraction of the dashed lines shown in panel (a). 
The black bar labelled "CI" indicates the position of the 
corresponding cluster as shown in Fig. [S] 

The lower panels of Fig. El show the histograms of the 
phase distribution in panel (b) for the Au film (c), the 
Au-film edge (d), the substrate (e), and the substrate 
edge (f). The full lines are Gaussian fits from Eq. Q 
except for panel (c) , where a double Gaussian of the form 

Wtot(A0) = OAu Wfilm(A0) + (1 - flAu) Wclus(A0) (6) 

with Gaussian contributions from the film Wfum and the 
adsorbed clusters Wcius was used. The width of the curve 
for the substrate of CTsubs — 0.18 rad is consistent with 
the averaged value obtained for a number of samples 
fsubs = 0.164 as discussed in Sec. IIII Al The difference 
may be attributed to poorer statistics of the relatively 
small sample size in Fig. [Sjand possibly to an increased 
number of pinned Au clusters near the islands. 

The substrate edge is expected to show growth pro- 
cess dependent steps which are yet too small to be re- 
solved in the noise of the density and vacuum fluctua- 
tions. Since a simple polynomial fit cannot account for 
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FIG. 6: (Color online) Details of scan 1 in Fig. Panel (a) 
shows the raw data of the phase shifts together with fitted 
values (dashed lines, see text), which average the Buctuations. 
Panel (b) shows the fluctuations Aij) obtained from the raw 
data after subtraction of the fits. Panels (c), (d), (e), and (f) 
show the binned data of the phase Buctuations form the Au 
Blm, the Au-film edge, the substrate, and the substrate edge, 
respectively. Solid lines are Bts from Eqs. Q and 0. 



this non-monotonous increase in thickness at the sub- 
strate edge, a broader distribution of the phase fluctua- 
tions is expected in that region. Indeed, the distribution 
for the substrate edge (Tsubs-cdgc = 0.24 [Fig. IHf)] is 30% 
broader than for the substrate plateau. Here a quadratic 
fit was appHcd to the substrate edge but Hnear or cubic 
fits give very similar results. 

The width of the phase distribution of the Au-film edge 
CAu-cdgc = 0.16 in panel (d) is close to the average value 
of the substrate, no smoothing of the surface corruga- 
tion is measurable. This result is consistent with the 
previously discussed interpretation (Sec. lIII^ . that the 
phase fluctuations of the substrate stem predominantly 
from intrinsic spatial density inhomogeneities. 

Finally, the bimodal distribution of the Au film in 
panel (c) of Fig. |H1 reflects the presence of adsorbed clus- 
ters on the substrate with an average height in this seg- 
ment resulting in an additional phase shift of A(/)cius ~ 
0.43 rad as indicated by the arrow. The width of the film 
part of the distribution aaim — 0.175 rad is comparable 
to that of the substrate and stems from the intrinsic spa- 
tial density inhomogeneities of the latter. For complete- 



ness we note the values of the cluster distribution width 
Ccius = 0.11 rad and the weighing factor aau — 0.695. 

C. Au-fllm thickness 

In principle the thickness of the Au film can be ob- 
tained from Eq. (Q. As discussed in Sec. Ill Bi thc param- 
eter of the mean inner potential Vq is not well defined as 
it appears to become thickness dependent for thin sam- 
ples. To be specific, Ichikawa and coworkers^*^ found for 
Au clusters on a Ti02 substrate with diameters smaller 
than 4 nm an increase of Vq of up to a factor of 3. In 
order to have a reference point for our system we have 
investigated the observed phase shifts of the cluster la- 
belled with "C2" in Fig.|21and compared the results with 
the shape of a cluster of the same size obtained in the 
MC simulations described in Sec. IV Bl 

Panels (a) and (b) of Fig. [7| shows the resulting phase 
shift profiles as obtained along the scans indicated in the 
inset. Dashed lines are linear fits. Panels (c) and (d) 
show the phase profiles of the cluster after subtraction 
of the linear contributions together with fitted ellipses 
(dashed lines). As a result of the strong fluctuations the 
fits are not unique. A width of 6cius = 2.3 ± 0.2 nm 
can be extracted while the phase shift is determined as 
A</)cius = 0.5 ± 0.02 rad. Comparing these values to the 
aspect ratio of ft-cius/^cius ~ 0.54 as obtained in Sec. IV Bl 
a height of /icius = 1-2 ± 0.1 nm is expected. Using Eq. 
n we obtain a mean inner potential of Vo — 57 ± 5 V 
consistent with the results from Ichikawa et alJ^ 

We have consequently confirmed the thickness depen- 
dence of the mean inner potential of Au clusters smaller 
than 4 nm in height. Unfortunately, when approxi- 
mating the apparent strong increase of the inner po- 
tential with decreasing thicknessiS in leading order as 
Vq ^ d^^ , the phase shift cx Vo d becomes indepen- 
dent of the thickness in leading order. From the phase 
shift induced by the Au film as measured in Fig. of 
A(/)fiini = 0fiim — '/'subs ~ 1 rad wc can thus only deter- 
mine an upper bound for the thickness of the Au film of 
c^fiim < 7 nm. This upper bound is obtained from Eq. 
using the bulk value'^^ Vau,gxp ~ 22 V, which is a lower 
bound for the inner potential. A more precise determi- 
nation requires the detailed knowledge of the sub- leading 
contributions to Vb(d). 

The intensity profiles in Figs. Qlb) and |S|^a) suggest 
that the film is thinner than the majority of the adsorbed 
clusters. This leads to an estimate of dfiim ~ 4 — 8 mono- 
layers. The observation discussed in Sec.^that clusters 
of size N < 300 (or dciustcr < 2 nm) are not stable enough 
to retain their structure on the substrate at room tem- 
perature and should consequently largely contribute to 
the film formation is consistent with these numbers. 



D. Heating effects 

The effects of heating on the shape of metal clusters 
has been studied early oniMiS In order to obtain infor- 
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FIG. 7: (Color online) Phase shift for two scans (a) and (h) 
through cluster "C2" in Fig.\^along the directions indicated 
in the inset. Dashed lines are linear fits to the background. 
Panels (c ) and (d) show the data of the cluster from panels (a ) 
and (b), respectively, after background subtraction. Dashed 
lines in (c) and (d) are elliptical fits (see text). 



mation about the stability of the islands we performed 
experiments heating the sample in situ. The sample was 
kept at 373 K for two hours. Figure IHl^a) displays the 
characteristics of a typical island before heating, whereas 
another island, which was not exposed to the electron 
beam during annealing, is shown in Fig. |SJb). Compar- 
ing these two images, it is obvious that as a consequence 
of the heating (i) the amount of Au forming the film is 
diminished, (ii) the number of particles is reduced and 
(lit) the sizes of the remaining particles have increased. 

The effects (ii) and (iii) are typical characteristics of 
Ostwald-ripening processes. The islands are metastable 
with an activation energy low enough that the applied 
heating of 100 K above room temperature is sufhcient to 
induce the transition on a time scale of a few hours. A 
more elaborate investigation of the dynamical behavior 
of the Au islands at elevated temperatures requires more 
experimental data and will be presented in a subsequent 
publication. 



V. MONTE-CARLO SIMULATIONS 

In order to obtain a more microscopical understand- 
ing of the formation of the Au islands — which cannot be 
resolved experimentally — we performed MC simulations 




FIG. 8: TEM images of two typical islands (a) before and (b) 
after in situ annealing at 373 K for 2 h. 

based on the Au-substrate interaction elaborated in Sec. 
mil and using the Gupta many-body effective potentials 
Eq. (O as described in Sec. IlIDI 

A. Cluster fusion and effective cluster potential 

As discussed in Sec. II V Al the observed Au film is likely 
to be composed of fused small Au clusters, which are 
present at the time of deposition but only unambiguously 
resolved in TEM for diameters larger than 1 nm (c.f. Fig. 

The fusion process of adsorbed clusters has been ob- 
served experimentally for Co clusters'* of diameter 8.5 nm 
as well as for clusters^ of Ausooo and and has been mod- 
eled for free clusters with MD methods using embedded 
atomSS and glue potentials^ for AU225 to Auasos- The 
objective here is to obtain an effective inter-cluster po- 
tential in the framework of the better suited many body 
Gupta potentials Eq. which we will use in Sec. IV CI 
to simulate the island formation. 

DFT results^^ show that the icosahedral configuration 
has a lower cohesive energy than the octahedral or cuboc- 
tahedral structures for the atomic closed shell cluster 
sizes N — 13, 55, 147. On the other hand, the Gupta 
potentials for Au exhibit a large number of local min- 
ima with amorphous structures and energies of AE^^ ~ 
0.01 eV lower than the icosahedral configuration.^^ Since 
AE55 - 0.01 eV < fcBTroom = 0.0256 eV, we use icosahe- 
dral configurations for the investigations presented here 
because different isomers are realized through thermal 
fluctuations at room temperature.^'* Note that only quan- 
titative details depend on the specific structure while our 
results are qualitatively quite general. 

Figure 1^1 shows the inter-cluster cohesive energy as a 
function of the distance of the center of mass points r 
of the two clusters. Panels (a) and (b) were obtained 
for two AU55 and two AU147 clusters, respectively. The 
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FIG. 9: (Color online) Fusion of Aum clusters for (a) N = 
55 and (b) N = 147. Gray and black circles are obtained 
in MC annealing runs at k^T = 0.0256 eV (T = 297 K) 
and ksT = 0.0001 eV (T = 1.16 K), respectively. The full 
lines are fits for the effective inter-cluster cohesive potentials 
Eq. Q). At low temperatures the fusion process comes to an 
early halt because the energy barrier involved in the atomic 
rearrangement cannot be overcome. 

gray points are obtained in free MC annealing runs at 
room temperature (fceT = 0.0256 eV) while the black 
circles were obtained at fceT 0.0001 eV (T = 1.16 
K). The scattering of the gray data points reflects the 
thermal activation of the clusters. At low temperatures 
the fusion process is halted at larger distances than at 
room temperature because the potential energy barrier 
involved in the reconstruction cannot be overcome. 

The full lines in Figs.|5Ia) and (b) are fits representing 
the effective inter-cluster cohesive energy and have the 
form 

^^AuTV = An {tanh [&A,(r - i?o,jv)] - 1} ■ (7) 

The parameters extracted are A55 — 2.702 eV, 655 = 
10.08 nm-\ i?o,55 = 0.524 nm and A147 = 4.941 eV, 
&147 = 7.149 nm~\ i?o,i47 = 0.718 nm. The fusion 
process stops at a temperature dependent minimal dis- 
tance i?min,Af(T') which Can be modeled in the poten- 
tial as a step function with a high repulsive value for 



r < Rmin,N{T), i.e., i^rcp = Acp8[Rmin,N{T) - r] with 
Aj-cp ^ Apf. The effective binding energy is then given 
by Ecs = EauN + Ercp. In the case of the investigated 
icosahedral structure the values are i?min,55('7room) ~ 0.4 
nm and -Rmin,i47(T'room) ~ 0.6 nm. Note that larger clus- 
ters are structurally more stable than smaller clusters 
and consequently the fusion process at U-qT = 0.0001 eV 
(T = 1.16 K) leaves the initial structure of the two AU147 
clusters more intact [Fig. Elb)] than in the case of two 
AU55 clusters [Fig. |5fa)]. 

B. Shape and adsorption energy of the adsorbed 
clusters 

For the interpretation of the holographic analysis of 
the Au film and the Au clusters it is useful to have in- 
formation about their shape. Moreover, from the simu- 
lations we obtain the effective attractive potential of the 
substrate exerted onto clusters of different size. 



1. Amorphous carbon substrate 

Since the exact shape of the substrate surface is not 
known (c.f. discussion in Sec. IIII|I we simulate the shape 
of clusters for flat surfaces with different relevant effec- 
tive mean substrate-adsorbate potentials derived in Sec. 
1111 Bl Figure [TUl shows a side view the shape of adsorbed 
clusters (substrate at the bottom) in the presence of a 
flat substrate with an attractive potential as given by 
Eq. O with a depth of yLj.o = 0-3 eV. The clusters of 
size (a) N = 55, (b) N = 147, and (c) N = 309 were 
prepared in icosahedral configurations close to the sur- 
face and subsequently freely evolved in a simulation run 
at room temperature (fceT = 0.0256 eV) for 10^ to 10^ 
Monte-Carlo steps per atom. 

As a result of the boundary condition imposed by 
the substrate the clusters reconstructed in the simula- 
tion process to hep (a) , fee (b) , and distorted icosahedral 
(c) structures. As expected, the larger N — 309 clus- 
ter is more stable than the smaller ones and does not 
undergo a structural transition within reasonable run 
times. The energy gain due to the Au-substrate inter- 
action is roughly the number of atoms at the surface 
times the potential depth I^lj.o- To be specific, we find 
A£;subs,55 = -6.31 eV, AS,ubs,i47 = -10.85 eV, and 
Ai?subs,309 = —12.02 eV for the configurations shown in 
Fig. ^1 as compared to the icosahedral structures. 

When briefly tempering the clusters with = 55 and 
N — 147 at k-QT ^ 0.7 eV and subsequent quenching we 
find that both hep and fee structures are realized. The 
energy differences of the hep configuration with respect to 
the bulk-ground-state fee structure is smaller than energy 
fluctuations due to boundary effects. This observation is 
consistent with the frequently observedS^ stacking faults 
in fee and hep bulk crystals. 

For completeness we have performed runs simulating a 
graphite substrate using Eq. Q with a potential depth 
of Vlj,o = 0-09 eV. We find that AU55 is amorphous at 
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FIG. 10: (Color online) Side view of the shape of adsorbed 
clusters of size (a) N = 55, (b) N = 147, and (c) N = 309. 
The initially icosahedral clusters are placed on a planar sub- 
strate in the x-y plane modeled with a generalized Lennard- 
Jones potential Eq.\^with a depth of the minimum of Vlj,o ~ 
—0.3 eV and simulated at room temperature (k^T = 0.0256 
eV). The clusters are reconstructed in hep (a), fee (b), and 
distorted icosahedral (c) structures as a consequence of the 
boundary condition imposed by the substrate. (The images 
were taken after quenching to T — 1.16 K in order to elimi- 
nate noise.) 

room temperature while AU147 retains its initial icosa- 
hedral structure (not shown). The melting point of 
Au55 in the absence of the substrate is close to room 
temperaturefSSiSLSiS242£ which accounts for its larger 
sensitivity to boundary effects. 



2. Gold film 

When placing a number of randomly distributed Au 
atoms on the homogeneous substrate as described by Eq. 
(O a hexagonally coordinated double layer is found to 
be the metastable structure attained (see also Sec. IV Dp 
for a large parameter range of 0.1 < li,j,o/eV < 0.6. 
The surface of such an Au double layer corresponds to 
an Au(lll) surface, albeit with 3.5% contracted nearest- 
neighbor bond-length due to surface effects 1^ In order 
to model the shape of the Au clusters that are found 
on the Au film we placed an initially icosahedral cluster 
N — 309 on the substrate surrounded by such an Au 
double layer formed of 291 atoms. The diameter of such 
a cluster of c? ~ 2.1 nm corresponds to that investigated 
in Sec. WU 

Figure ITTI shows the resulting shape of the cluster from 
the side in panel (a) and from the top in panel (b) . Since 
the Au-double-layer-cluster interaction is much stronger 
than the substrate-cluster interaction [c.f. Fig. llOf c')]. the 
cluster completely reconstructs to an fee lattice to match 
the structure of the film. 



C. Island formation 

In order to model the formation of the islands on the a- 
C substrate we simulate the dynamics of clusters rather 
than individual atoms, which allows us to model suffi- 



FIG. 11: (Color online) Side (a) and top (b) view of an AU309 
cluster surrounded by an Au double layer bound to a substrate 
as modeled by Eq. |^ with Vl,t,o = 0.3 eV at room tempera- 
ture (kuT = 0.0256 eV). The AU309 was initially icosahedral 
but reconstructed to a fee lattice under the inHuence of the 
Au(lll) double layer. (The images were taken after quench- 
ing to T = 1.16 K in order to eliminate noise.) 



ciently large systems. Moreover, we are able to reproduce 
experimental result a^i^i^ for the diffusion and agglomer- 
ation of size selected clusters published previously. 

In Sec. IIV CI we have shown from the experimental 
holographic data that the Au film in the islands is less 
than 7 nm. In Sec. IV Bi the MC simulations have shown 
that small adsorbed clusters of Auat for = 55 and 
N = 147 do not retain their structure at room tempera- 
ture and can easily fuse. Smaller clusters are more likely 
to be trapped by pinning centers on the corrugated a- 
C film as discussed in Sees. IIIl'Xl and ITlI Bl while larger 
clusters are more stable and retain their structure instead 
of fusing with their environment as indicated in Sees. IV A1 
and lVBI These findings suggest that the Au films of the 
islands are formed by diffusing and coalescing clusters on 
the substrate in the size range around 50 < < 300 
leading to a film thickness of 4 to 8 monolayers. 

For simplicity we show here simulations of a sys- 
tem with an effective cluster cohesive potential of the 
form given by Eq. lO with parameters as obtained for 
A^ = 147. Correspondingly, the cluster-substrate attrac- 
tive potential is given by Eq. (jSJ, where we use a value 
of Vlj,o = 11 eV as obtained from the quantitative anal- 
ysis in Sec. IV Bl Results obtained for parameters with 
A^ = 55 are qualitatively equivalent (not shown) . 

Figure IT^ shows a series of snapshotsSi from a simula- 
tion run of 200 initially randomly placed effective AU147 
clusters. Panel (a) shows that after 3000 MC steps per 
cluster small islands and chains of clusters have formed 
very similar to those observed on short time scales for 
monodisperse 4.8 nm diameter Au clusters on an a-C 
substrate^ Agsooo clusters on graphite;^ and 8.5 nm di- 
ameter Co clusters on microgrid substratesi^ It is impor- 
tant to note that the time scales for the diffusion and fu- 
sion processes depend on the temperature, the substrate 
properties, and the cluster sizes so that these metastable 
structures can be observed only in a system specific time 
window. 

At room temperature the AU147 clusters can be consid- 
ered as to be fused almost instantaneously compared with 
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FIG. 12: (Color online) Simulation of 200 AU147 clusters in- 
teracting as described by Eq. on a substrate modeled by 
Eq. 0) with Vl],o = 11 eV after (a) 3000, (b) 4 x lO"*, (c) 
4 X 10^, and (d) 4.3 x 10^ MC steps per cluster. The ra- 
dius of the dots representing the clusters has been chosen as 
,i47(rroom) = 0.6 nm, i.e., toucMng clusters are fused 
as shown in Fis. IPKb). The side views in panels (a) and (b) 
show that all clusters are located in a plane touching the 
substratei— 



the diffusion time scales as shown in Fig. |5^b). Corre- 
spondingly the radius of the dots representing the AU147 
clusters in Fig. 1121 has been chosen as -Rmin, 147(7100111) = 
0.6 nm. Larger clusters fuse on much longer time scales 
or at higher temperatures From the MC simula- 
tions a diffusion coefficient cannot be determined quan- 
titatively because of the lack of a time constant in the 
MC procedure. 

Panels (b), (c) and (d) of Fig.lT^shows the system after 
4x10^^, 4x 10^, and 4. 3x 10"^ MC steps per cluster, respec- 
tively. Clusters and islands of clusters coalesce once they 
get close enough gaining cohesive energy. Structures with 
short border lines are favored over extended linear ones 
for the same reason. The (expected) tendency to form 
larger islands by diffusion and fusion is clearly visible. 
The present simulation of effective clusters underesti- 
mates the migration speed as well as the reorganization^^ 
of the islands to more circular structures because the dif- 
fusion of the individual atoms on the surface of the clus- 
ters is not accounted for. For larger clusters the latter 
is the dominant cluster migration process^'^^ referred to 
as surface diffusion. Including this effects requires the 
investigation of much smaller systems such as shown in 



Sec.rVDl 

In conclusion the cluster diffusion and fusion as mod- 
eled by the effective clusters and cluster interactions cor- 
rectly accounts for the experimentally observe d^i^'^ pat- 
tern formation [Fig. ll2r a')]. The fused clusters [Fig.|5fb)] 
form small islands [Fig. ll2r d)] with a thickness of roughly 
the cluster diameter, which corresponds in the case of 
Aui47 to 5 monolayers [Fig. [rnT b)] . The expected con- 
tinued migration by surface diffusioni*^^ and fusion of the 
small islands in Fig. I12f d) is consistent with the experi- 
mentally observed island formation shown in Fig. ^b) . 



D. Cluster drag 

All initially randomly distributed larger clusters are 
incorporated into the Au films in the islands as shown 
in Figs, n^b) andlH^a). It is not possible to simulate 
this effect directly because of the size limitations of the 
method. It is possible though to place randomly dis- 
tributed Au atoms on a substrate modeled by the gener- 
alized Lennard Jones potential Eq. where T^lj,o = 0.3 
eV, together with a N = 147 and a = 55 Au cluster at 
room temperature {k^T = 0.0256 eV). The thus modeled 
system is significantly smaller than those observed in the 
experiments but the dynamical behavior is likely to be 
similar, albeit with rescaled diffusion constants and time 
scales. 

Snapshots taken at three different run times are shown 
in Fig. El Panel (a) shows the initially random distribu- 
tion of 298 atoms as well as the AU55 and AU147 clusters 
after 93 MC steps per atom. Panel (b) shows the system 
after 6 x 10^ MC steps per atom. The percolating regions 
coalesce. Panel (c) shows a side view of the same state of 
the system as in (b) and illustrates that the atoms have 
formed a Au double layer. This double layer is stable 
at room temperature in a substrate parameter range of 
0.1 < yLj.o/eV < 0.6. 

Panel (d) of Fig. Ill shows the system after 1.18 x 10^ 
MC steps per atom. The coalescing Au double layer has 
contracted thus minimizing it cohesive energy. In the 
cohesion process the less mobile AU55 and AU147 clus- 
ters have been pulled closer together. They are located 
closer to the border of the Au film consistent with the ex- 
periments indicating a slightly larger probability to find 
clusters near the boundaries of the islands than in their 
center. The separated island in the lower right hand cor- 
ner of Fig. El is less mobile and fuses with the larger 
island only after comparably long run times (not shown) . 

Panel (e) of Fig. El shows a side view of the system 
after 1.18 x 10® MC steps per atom and illustrates the 
fee restructuring of the initially icosahedral clusters in 
the presence of the Au double layer film. 



VI. CONCLUSIONS 

We presented experimental results showing the forma- 
tion of Au islands on amorphous carbon substrates after 
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FIG. 13: (Color online) Time evolution of a system of 298 Au 
atoms placed randomly on the substrate [Eq. {^], one AU55 
and one AU147 cluster. Panel (a): 93 MC steps; panels (b) 
and (c): 6 x 10* MC steps; panels (d) and (e): 1.18 x 10'^ 
MC steps per atom. Panels (c) shows the formation of the 
double layer film, panel (e) the restructuring of the initially 
icosahedral clusters. (The images were taken after quenching 
to T — 1.16 K in order to eleminate noise.) 

laser deposition of non-size-selected clusters and subse- 
quent aging of the samples in the absence of inert con- 
dition for three to four months. We characterized both 
the substrate and the Au films with electron transmis- 
sion holography. A model potential for the substrate- Au 
interaction was derived with the help of density func- 
tional calculations. Subsequently a number of Monte- 
Carlo simulations were carried out describing the island 
formation process and the shape of the adsorbed clusters. 

The following qualitative and quantitative results have 
been incurred from the investigations. 

• The observed islands have a typical size oi A ^ 
25 X 70 nm2. 



• The islands consist of an Au film of a few mono- 
layers [probably 4 to 8 but definitely less than 27 
Au(lll) layers] thickness and adsorbed larger clus- 
ters. 

• The islands are metastable and coalesce at elevated 
temperatures of T ~ 400 K. 

• The islands are formed by diffusion and coalescence 
of clusters of size 50 < N < 300. 

• Larger clusters, which normally are less mobile, are 
dragged along by the percolating Au film formed by 
the smaller clusters. 

• Smaller clusters are pinned at deeper substrate in- 
trusions as indicated by the broadening of the phase 
distribution in the holographic images of the ab- 
lated samples. 

• A number of observations suggest that the amor- 
phous carbon substrate surface is much less corru- 
gated than anticipated. Instead, the substrate ap- 
pears to exhibit significant internal spatial density 
fluctuations. 

• The Au-substrate interaction can be modeled by 
a generalized Lennard- Jones potential for graphite 
[Eq. |@J] and for amorphous carbon [Eq. ©]• 

• An effective Au cluster-cluster interaction has been 
determined [Eq. Q]. 

The extensive investigation presented in this work con- 
tribute to a better understanding of some processes tak- 
ing place in the so far little investigated field of the dy- 
namic behavior of cluster arrays deposited on amorphous 
surfaces. 
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